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Abstract 

We study non-Abelian solitons of the Bogomol'nyi type in M = 2 (d = 
2 + 1) supersymmetric Chern-Simons (CS) and Yang-Mills (YM) theory with 
a generic gauge group. In CS theory, we find topological, non-topological and 
semi-local (non-)topological vortices of non-Abelian kinds in unbroken, broken 
and partially broken vacua. We calculate the number of zero-modes using an 
index theorem and then we apply the moduli matrix formalism to realize 
the moduli parameters. For the topological solitons we exhaust all the moduli 
while we study several examples of the non-topological and semi-local solitons. 
We find that the zero-modes of the topological solitons are governed by the 
moduli matrix Hq only and those of the non-topological solitons are governed 
by both Hq and the gauge invariant field f2. We prove local uniqueness of the 
master equation in the YM case and finally, compare all results between the 
CS and YM theories. 
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1 Introduction 

Chern-Simons (CS) theory is an alternative gauge theory in odd space-time dimensions (we only 
consider d — 2 + 1 here), which is very different from Yang- Mills (YM) and Maxwell gauge 
theories. The CS gauge fields in pure CS theory are non-dynamical, however, they become very 
interesting upon coupling with other fields. For instance, when considering a CS kinetic term 
together with a YM kinetic term, the gauge bosons acquire a topological mass in the absence 
of the Higgs mechanism pQ. The CS gauge theories are interesting both due to their theoretical 
beauty and also for their experimental applications, such as the fractional quantum Hall effect 
The CS term can be induced by quantum radiative corrections, even if the original Lagrangian 
does not include a bare CS term [3] , such as the d = 2 + 1 dimensional quantum electrodynamics 
(QED). 

One of the important features of both CS and YM theories coupled with scalar fields is the 
existence of solitonic solutions to the classical equations of motion. Thousands of works are 
made on the study of these solitons. For the CS solitons, early works with/without the Maxwell 
term were studied in Refs. [U [5] . In the Abelian CS theory coupled with a complex scalar field, 
solitons of the self-dual type, viz. solutions satisfying the Bogomol'nyi energy bound and in turn 
requiring a special sixth-order scalar potential were found in Refs. [HJ [7] . The supersymmetric 
version which corresponds to the self-dual model was obtained in Ref. [8] and the study on zero- 
modes of the solitons was developed in Ref. [9]. Furthermore, due to the presence of both the CS 
vacuum and the Higgs vacuum, domain walls and non-topological vortices were found in addition 
to the familiar topological vortices in the Abelian CS-Higgs (CSH) model [9]. The non-topological 
soliton is a typical soliton in CS theory which in fact does not exist in the Maxwell-Higgs (or 
Yang-Mills-Higgs (YMH)) theories. Similar solitons of topological/non-topological kinds were 
found also in the Maxwell-CS-Higgs (MCSH) theory [101 [Til EE] • Furthermore, a semi-local 
extension of this type of solitons, viz. in MCSH coupled with two Higgs fields was found both in 
the topological as well as the non-topological case in Ref. [T5] . An extension to the non- Abelian 
(SU(N)) CS theories was also made and non-topological solitons with a global charge were found 
in Ref. [H]. There are also lots of interesting works on the solitons in the non-relativistic CS 
theories. Many relevant references on both the relativistic and non-relativistic models can be 
found in the excellent reviews [T5| [To] . The mentioned works on the CS solitons were established 
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in the early '90s and we shall refer to all these as solitons of the Abelian kind 

The second era of the CS solitons recently started when a new type of the toplogical CS 
soliton of the non- Abelian (NA) kincfl was found in the U(N) CS theory coupled with N Higgs 



fields in the fundamental representation [21]. The NA CS vortex found so far is quite similar to 
the NA YM vortex [TTI [TBI [T9l 120] . It lives in the color-flavor locked, broken vacuum (i.e. the 
Higgs vacuum) and has internal orientational zero- modes, CP^ -1 . These internal orientational 
modes are related to the color and flavor degrees of freedom, which make the solitons being truly 
of the NA kind. Soon after its discovery, the NA semi-local vortices were found in the U(N) CS 
theory with N{ > N fundamental Higgs fields [221 123] . The dynamics of the NA vortices was 
studied in Ref. [21] and the NA vortices with a globally conserved charge was found in Ref. [2B] . 
The topological solitons of both local, semi-local and fractional kind in the NA CS theory with 
an arbitrary gauge group of the form 17(1) times a simple group have also been studied very 
recently in Refs. [261 [27] . Thus, the study on the CS solitons is being revived now. 

The NA CS solitons found so far are all topological solitons in the broken vacuum (the Higgs 
vacuum). This is because their discovery [2T] was inspired by the NA YM vortex which is indeed 
a topological soliton living in the broken (Higgs) vacuum. However, we know that there are also 
non-top ological solitons in the Abelian CS models [TJ- So it is natural to ask ourselves if non- 
topological or other solitons of the NA kind exist or not, in NA CS theories. With this question 
in mind, we will be concerned in this paper with vortices of all possible non- Abelian kinds in 
the NA CS gauge theories with the gauge group G being general; G = U(l) x G' with G' is an 
arbitrary simple group. As if it was not broad enough a spectrum, we will make every statement 
throughout the paper be independent of the gauge group. Beyond our naive expectation, we will 
find semi-local NA vortices of both topological and non-topological types in the models where 
the number of flavors is less than or equal to the number of colors. 

We will investigate not only single solitons but also multiple solitons, viz. we will focus on 
the zero-modes of the Bogomol'nyi solutions. To this end, we will count the number of zero- 
modes by using an index theorem technique. Then we will extend the moduli matrix method 
for explicit realization of the counted moduli. The moduli matrix formalism was first introduced 
in the studies of domain wall systems [281 12H] but later developed to encompass vortex systems 



In CS theory as well as YM theory, the vortex in SU(N) theories is basically of the Abelian kind due to the 
first homotopy group characterizing the vortices being -ki{SU{N)/'Lm) = and hence these vortices do not 
carry non- Abelian orientational moduli, as opposed to those we discuss next. 

2 The non- Abelian vortices were first found in the supersymmetric YM theories in four dimensions [T7| ITS] . 

There are many works on the NA vortices in YM theories and we shall refer to the literature for the many results 
which however are nicely summarized in several review articles |191 120) . 
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[30] . which has proven an invaluable tool in soliton studies [3TJ [32j (33J [20]. The moduli matrix 
formalism has already been introduced to the CS models in Refs. [261 I2Z] for the topological 
solitons. Here we will try to further develop it to cover all the solitonic solutions in the system, 
including the non-topological solitons in the unbroken and partially broken vacua. It turns out 
that the moduli matrix formalism, which has proven itself incredibly useful for topological solitons 
in YM theories, provides us with a common framework for studying all kinds of solitons not only 
in YM but also in CS theories and especially in all vacua of the systems at hand. In order to 
find similarities and differences among the solitons in YM and CS theories, we will show all the 
results parallelly in the two cases. Most of the statements concerning the YM vortex (up to 
Sec. H]) are just a review of the known results, however written in a more transparent way with 
respect to the gauge group. A new result - and one of the main results of this paper - for the 
YM vortex is on the uniqueness problem of the so-called master equation (Sec. |5]). 

As we have already mentioned, we will pay special attention to the moduli space of solutions 
(or simply moduli space). It is one of the most important properties of solitons - of topological 
or non-topological kind. The moduli space has a certain dimension characterizing the soliton 
and the gauge theory (or string theory configuration) at hand and some local coordinates on 
it which can be identified with the bosonic zero-modes of the system describing the soliton. 
The number of bosonic zero-modes gives the aforementioned (complex) dimension of the moduli 
space. The moduli space per se is a classical statement which eventually, far in the infrared, will 
be modified by quantum corrections. However, this information usually is available only through 
the knowledge about the classical moduli space supplemented by calculations of a corresponding 
sigma model describing the low-energy effective theory of the soliton under consideration. 

A further powerful technique often used in soliton physics is the moduli space approximation 
or the so-called geodesic approximation due to a seminal paper [M] which was first applied 
to monopole scattering and thereafter to a vast number of other soliton configurations in the 
literature. 

Index theorems have proven to be very powerful tools in physics. Needless to mention is the 
Witten's index indicating if a theory at hand can have its supersymmetry broken spontaneously 
or not [35] . There is also the Hopf index theorem in sigma models giving the number of vacua in a 
class theories [36]. An index closely related to what we are interested in here, is the Atiyah-Singer 
[37] index theorem, which has been used to count the physical parameters describing instantons 
after the configuration space has been compactified [38J. Finally, there is a generalization of 
this index to non-compact manifolds, namely the Callias index theorem which counts the zero- 
modes of a Dirac operator minus the zero- modes of the corresponding adjoint operator [39] . 
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This is a very useful technique in soliton physics and especially in BPS systems where the 
fluctuations of the BPS equations can easily be written as a Dirac operator acting on a (vector) 
space of fluctuations. To mention a few cases in the literature, this technique has been applied 
to monopoles [10] and Abelian-Higgs vortices [H]; to Abelian CS vortices [9]; to Maxwell-CS 
vortices [TT| ; to non- Abelian U(N) vortices [I7J and to domain walls of Abelian kind [12] and 
non- Abelian kind |43j . 

The organization of the paper will be as follows. In Sec. [2] we set up our theories, the YM- 
Higgs model and the NA CS-Higgs model. Here we will review the basic properties including the 
BPS-equations of both theories. In Sec. El we will derive a generalized formula for the Callias- 
type index of a certain class of BPS systems including both our models under consideration. In 
Sec. H] we will review and develop the moduli matrix formalism to realize the moduli parameters 
and also explain the new types of solitons, namely the NA non-topological solitons. In Sec. [5] 
we will confront the long-standing problem of the uniqueness of the master equations. We find 
a relation between the variation of the master equations and the vanishing theorems studied in 
Sec. El We conclude the paper with a discussion and outlook in Sec. El Sees. 13.21 and 14.11 are 
reviews of known results, so the reader who is familiar with the NA YM vortices can skip them. 

2 The model and notation 

We begin with the M = 2 supersymmetric Yang-Mills-Chern-Simons (YMCS) theory coupled 
with Higgs fields in d = 2 + 1 dimensions. In order to make the following arguments applicable to 
a wide class of gauge theories, we will not specify the gauge group unless we make some explicit 
examples. Indeed, we will take the gauge group to be on the form G = {U(l) x G') /Z no , where 
G' is always a simple group. When we will make some examples we will use the gauge group G 
with G' = SU(N), SO(N), USp(N). In this case we can choose either N = 2M or N = 2M + 1 
for SO(N) whereas N = 2M for USp(N). Here Z no is the center of G', see Table EJ 



G' 


SU(N) SO{2M) USp(2M) SO{2M + 1) 


n 


N 2 2 1 



Table 1: Various values of the greatest common divisor (gcd) of the Abelian charges of all the G' 
invariants which have a non-zero VEV at infinity (in the Higgs phase). 
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We will use the standard convention for the Hermitian generators 



1 



Tr(t a t?) = -5 a)3 , t' 



for a,P = 0,1,2, 



2 V2iV 
dim(G'), while the index of G' will be denoted by a, b = 1, 2, 



(2.1; 



dim(G'). 



The G' vector multiplet contains the gauge fields and a real adjoint scalar field {A^, a }. The 
U(l) vector multiplet also has the corresponding gauge fields as well as a real singlet scalar field 



{A , 0°}. We consider N f Higgs fields H J T 



1, 2, . . . , N; A = 1,2, ... , Nf) in the fundamental 



representation of G' with uniform U(l) charge. In the following we will use a matrix notation 
where H is an N x Nf dimensional matrix, so that the gauge symmetry acts from the left-hand 
side and the flavor symmetry acts from the right-hand side. We will use the following compact 
notation for the construction of g and g' algebra valued fields, respectively, as follows 



dimG' 

E ■ 

a=0 



dimG" 

E ' 

a=l 



(2.2) 



We are now ready to write down the Lagrangian. Since all the fermions do not play an 
important role in the following argument for the solitons under consideration, we show only the 
bosonic part 



C 



YMCSH 



4e 2 
+ Tr 

£_ 

2 



2e 



8vr 



Wmcsh — — \ Tr 
Our conventions are 



V^H (V»H) ] 



27r' 



Wmcsh 
t a 



e 
~2 



Tr 



(hh* 



K 

2^' 



(2.3) 
(2.4) 



d^A u — d u A^ + i [Ap, A L 



iAjH 



V lt <f> = d li <f> + i[A^,<l>] . (2.5) 



The first line is the kinetic term of the G 1 vector multiplet and the non-Abelian CS term with 
the gauge coupling constant g and the CS coupling constant fi which has to be an integer in the 
quantum theory in order to preserve gauge invariance (up to large gauge transformations). The 
second line is the contribution from the Abelian vector multiplet with e being the gauge coupling 
constant and k being the CS coupling constant which can take on any (real) value. The potential 
consists of two terms. The first term is the D-term of G' vector multiplet while the second term 
is that of the U(l) vector multiplet. 
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Although we have turned off the terms including fermions, one can still smell the super- 
symmetric nature in the special relation between the gauge and scalar coupling constants. The 
so-called Fayet-Iliopoulos (FI) parameter £ has been chosen to be positive which ensures stable 
supersymmetric (SUSY) vacua. In order to simplify the notation we define also the parameter 

| > . (2.6) 

Note that the model has SU(N{) flavor symmetry. 

If we completely discard the Higgs fields H (i.e. setting H = 0) of the full Lagrangian ( 12. 4p . 
the FI term in Eq. ( 12.4ft can be absorbed by a constant shift in 0°. Hence, the vacuum = is 
in the symmetric phase where no symmetries are broken. Although the Higgs mechanism is not 
at work, the vector multiplet acquires a topological mass which is given by 

Ke 2 [xg 2 , , 

Here m K is the mass of the Abelian gauge fields A ^ as well as the real scalar field 0° while m M is 
the mass of G' gauge fields as well as the real adjoint fields a . 

In the following we will study the above described model in two limits, namely k, fi — > which 
reduces the model to the Yang-Mills-Higgs (YMH) theory without the CS interactions and e, (7 — 
00 which reduces it to the "pure" Chern-Simons-Higgs (CSH) theory without dynamical gauge 
fields. We will leave the intermediate case of both kinetic terms at finite coupling, viz. YMCS- 
Higgs theory for a companion paper 



2.1 Yang-Mills-Higgs theory 



Taking the limit eliminates the CS kinetic term for the gauge fields and decouples the 

adioint scalar field from the D-terms. Hence the bosonic Lagrangian density can now be written 



C 



1 



1 



YMH 



+ Tr 



V ra 2 + 27 ( ^ )2 -i ra2 '2e 2 

92 {Tr [HHH a ] } 2 - - {Tr [(Hrf - v 2 l N ) t°] }" . (2.8) 



2 1 L JJ 2 

In the case of G' = SU (N) there is a unique Higgs vacuum where G is completely broken 

H = vl N , = 0. (2.9) 



3 This model can be trivially embedded in 3 + 1 dimensions where the vortex solutions describe strings instead 
of particles in which case the supcrsymmetry allowed by the given potential is J\f = 1 in 3 + 1 dimensions. Either 
point of view is consistent with our discussion in the following, i.e. formally there will be no difference. 
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The Higgs phase is in the color-flavor locking phase where the diagonal global symmetry is 
preserved 



U(N) x SU(N { ) -> SU(N) c+{ . (2.10) 

The vacuum is gapped and there are no massless modes, while the mass spectrum is given by 

m e = ve , m g = vg . (2-H) 

Thanks to supersymmetry, the vector multiplets and the chiral multiplets have the same masses, 
viz. the Abelian gauge fields the real scalar field cjP and the real part of the trace part of 
H all have the same mass m e . On the other hand, the SU(N) gauge fields A®, the real adjoint 
fields <p a and the real part of the traceless part of H all have the mass m g . 

For G' = SO(N), USp(N) the vacuum structure is much more complicated, see Ref. [45] . We 
will however, only consider the same vacuum ( 12.91) for the reason that it preserves the maximal 
global color-flavor locking symmetry 

G x SU(N { ) -> G' c+f . (2.12) 

The masses of the vector multiplets remain the same as in the G' = SU(N) case. Those for the 
trace part and some of the traceless part of H are unchanged but the rest of the fields contained 
in H become massless (i.e. due to flat directions) since the number of gauge fields is not sufficient 
for the Higgs mechanism to eat all the massless fields. 

Let us consider topological solitons in this model. We are interested in the vortex which is 
usually called the non- Abelian vortex. It is a natural extension of the Nielsen- Olesen vortex of the 
Abelian-Higgs model. Performing a Bogomol'nyi trick on the Hamiltonian under the assumption 
that the configurations are static, i.e. the energy in the C-plane 

T= J c (± (*£ - 9 2 Tr [HHh«]f + ± « - e 2 Tr [( H rf - vH N ) t ]) 2 (2.13) 

+ ^(W) 2 + ^i0°) 2 + Tr [4 \VH\ 2 - v 2 F 12 - ie«0, {(V 3 H) H^]\ \ , 

we can read off the BPS tension which is the lower bound 

T B ps = -v 2 / Tr [F 12 ] = 2ttv 2 Nu = 2<z/ > . (2.14) 

Here v is the Abelian winding number determined as 

" = -Af /Tr[F 12 ] = - , fceZ >0 , (2.15) 



where n is the center of the gauge group Z no (see Table Finally, from Eq. ( 12. 13f) we have 
the BPS equations (<p a = 0) 

VH = 0, F£ = g 2 a Tr [(HH* - v 2 l N ) t a ] , (2.16) 

where g a stands for the gauge coupling constants (go = e and g a = g) and a is not summed over. 
The non-Abelian part of the second equation is expressed in matrix notation as 

F 12 =g 2 {HH^) G , , (2.17) 

where we have introduced a bracket for a projection operation 

(X) G , = Tr [Xt a ]t a , (2.18) 

for an arbitrary N x N matrix X. 

In the special case of e = g we can simplify the second BPS equations to the following 
combined equation 

F 12 = e 2 (HH^ -v 2 l N ) G , (2.19) 

where we have defined 



(X) G ee Tr [Xe] e = ^Tr[X]l w + (X) G , . (2.20) 



Note that (X) G = X/2 holds HXeG. 

The above special case of e = g is not only simple but also has the advantage that solutions 
to the Abelian BPS equation automatically solve the non-Abelian BPS equations. Hence, with- 
out solving the non-Abelian equations, we already have a known class of solutions which are 
essentially Abelian embeddings. All the minimally winding solutions are indeed in this class. Of 
course, we have to solve the non-Abelian system for the higher winding vortices which are out of 
the class. 

In particular, we can write the explicit form of (X) G , with G' = SU(N), SO(N) and USp(N) 

as 

G^SUiN), 

G \ i (X - JtA T j) , G> = SO(N), USp(N) . 



4 no is the greatest common divisor (gcd) of the Abelian charges of all the G invariants which have a non-zero 
VEV at infinity in the Higgs phase, see Ref. [46] . 
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For U G G' = SO(N), USp(N) we will use the invariant rank- two tensor J defined as 

U T JU = J , J f = eJ , JKJ = 1 N , 



(2.22) 



with e = +1 for SO(N) and e = — 1 for USp(N). Throughout the paper we will adapt the basis 
in which J is 




J 



for SO(2M) (e = +1) and USp{2M) (e = -1) while for SO(2M + 1) instead 



(2.23) 





( 


1a/ 


o\ 


J = 


1a/ 










1° 





V 



Let us give some examples of minimal winding solutions in the case of G' 
and f7Sp(4) 



5t/(4) 
SO(A) t USp(A) 




H = diag 



) 



= diag (i^NO^ o,0,0) 
diag 



ANO r/ANO 
fe=l ) - n fc=l ) 



V, V 



) 

diag (F^i^i, 0,0) 



(2.24) 
SU{4),SO{4) 

(2.25) 
(2.26) 



where (if^ NO , -F^ ) stands for fc G Z> winding coaxial-vortex solution in the Abelian theory. 
Generic non-Abelian solutions in the 577(4) and USp(4) models can be generated by acting with 
the global symmetry of the system on the above solutions. There exist further solutions which 
belong to the different classes in the SO (A) cases. 

The most significant feature of the non-Abelian vortex is the presence of the so-called orienta- 
tional zero-modes associated with the spontaneous breaking of the color-flavor locked symmetry. 
In the above example of G' = S77(4), one can easily see that S77(4) c +f obeyed in the vacuum 
(r — > oo), is in fact broken to C/(3) c +f at the center of the vortex where 77^° — > 0. Thus the 
Nambu-Goldstone zero-modes of the single non-Abelian vortex in the G' = SU (4) theory are 



M k=1 

JVl SU{A) 



C x CP 3 



(2.27) 



Here the first factor, C, corresponds to the position zero-modes. Similarly, for the case of G' 
USp{A), the vortex breaks USp(4) c+ t down to Z7(2) c+ f . Hence, in this case the zero-modes are 

USp{4) 



C x 



U(2) 



(2.28) 
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while for G' = SO (4) there exists a topological charge in addition to the vortex number, which 
is called the Z 2 charge [47j. The moduli space is then split into two copies of the space (HI 147] 
given by the breaking of SO(4) c+ i down to U(2) c+ { as 



•^SO(4) 



C x 



50(4) 



17(2) 



U 



50(4) 



17(2) 



(2.29) 



In the case of G' = SO(N), with N > 4 the zero- modes are somewhat more elaborated, so we 
will not explain the details here, which however can be found in Ref. [47] . 



2.2 Chern-Simons-Higgs theory 

The second limit under consideration is e, g — > oo, i.e. the strong gauge coupling limit where the 
kinetic terms of the vector multiplets vanish. Then the original Lagrangian density (12.41) reduces 
to the CSH model (see Ref. 



+ Tr [V^H - Vcsu , (2.30) 

^csh = Tr [0 2 Mt] . (2.31) 

The sixth-order scalar potential Vqsh is given by eliminating the adjoint scalar fields <p a from 
Eq. (EZD by 

K r 



-Tr [(M* - v 2 l N ) t a ] , a = 0, 1, 2, . . . , dim(G') , (2.32) 



with kq = k and K a — for a — 1, 2, . . . , dim(G'), and a is not summed over on the right-hand 
side. 

There exist three types of vacua in this model: 

• Symmetric phase: One vacuum is in the completely symmetric phase (the CS phase) where 
(H) = and no symmetries are broken. The vector multiplets are decoupled, so the Higgs 
fields are the only dynamical degrees of freedom with the mass 

m H = m = . (2.33) 

• Asymmetric phase: There is also a vacuum in the Higgs phase (H) = uljy, where the gauge 
symmetry is completely broken. The mass of Higgs fields is the same as that of gauge fields 
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due to supersymmetry. The Abelian gauge fields A° u and the real part of the trace part of 
H have the same me 
H have the mass m, 



H have the same mass m KOD , while the G' gauge fields A" the real part of traceless part of 



"/JOO- 

ml Attv 2 m 2 Airv 2 . A . 

m KOO = — = , m MOO = — = . (2.34) 

The Higgs vacuum {H) = vIn is unique if G' = SU(N) but there exist flat directions in 
the cases of G' = SO(N), USp(N) as we have seen in the previous section. 

• Intermediate phases: In between there is a variety of partially broken phases. 

In the model with G' = SU(N) and Nf = N, the vacua are labeled by an integer m = 
0,1, 2,..., N as follows 



H {m) = diagO, v, ■ ■ ■ , v, 0, ■ ■ • , 0) . (2.35) 




The vacuum m = N is the full Higgs vacuum where the symmetry U(N) C x SU(N){ is broken 
to SU(N) c+ {. The vacuum m = is the unbroken phase where no symmetries are broken, while 
in the case ofl<m<iV — lwe have the intermediate vacua. In this case the symmetry is 
broken down to U(N — m) c x S[U (rra) c +f x U (N — m)f]. Note that the vacua with m = 0,N are 
unique but the rest are not. The vacuum manifold is a complex Grassmannian manifold given 
by 3U(N) t /3[U(m)^ f X U(N - m) f ] ~ Gr N>m . 

The variety of the vacua results in a variety of vortices which can be either topological or 
non-top ological. In order to derive the BPS equations, we again perform a Bogomol'nyi trick on 
the tension of the vortex 

T= f Tr \\(V -i(f))H\ 2 + A\VH\ 2 -v 2 F 12 -ie ij d i {(V j H)H^}] , (2.36) 

where is given by Eq. (I2.32p . The tension for the BPS saturated vortices is given by 

T BPS = -v 2 / Tr[F 12 ] = 2ttNv 2 u > . (2.37) 
Jc 

where v is the U(l) winding number 

for topological solitons 

for non-topological solitons , (2.38) 
^ + for topological and non-topological solitons 

where k, k' are integers and a is a real number. In the case of Abelian non-topological solitons, 
the integer part k is limited as a > k + 2, see Ref. [19]. For non- Abelian cases, no conditions 



2vrV2iV Jc 
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for a and k are in general known yet. But only in the special case of U(N) gauge group where 
the coupling constants are equal k = fi and for the diagonal solutions, we can easily extend the 
Abelian results as a« > ki + 2. Here ccj (a = and ki (k = kj) are associated with the 

non-topological and topological flux of the i-th U(l) subgroup of U(1) N C U(N). 

Let us furthermore define the magnetic flux as follows 

$ = 2nNv . (2.39) 

The BPS equations in this case read 

VH = , V H = i(pH , (2.40) 

which however due to the presence of electric charge density in the vortex have to be accompanied 
by the Gauss law (a not summed over) 

((V H)Ht-H(V H)^t a ] . (2.41) 
Combining the equations (I2.40p and (12.41 p we obtain the following system 

A-TT 

VH = , F? 2 = — Tr [{0, HH ] ) t a ] . (2.42) 

This system is valid for any simple group G' . In matrix notation, we can explicitly write down 
the second equation as 

= [HH< - vH N ] Tr [HH<] 1 N + ^Tr [(HH% HH>] 1 N , (2.43) 

where we have used the following identity 



12 ~ K 



4-bl +a 



(X (X)) G , = Tr [XtH a ] Tr [Xt b ] t a 

= Tr [Xt a t b ] Tr [Xt b ] t a + i/ bac Tr [Xt c ] Tr t° 
= ((X)X) G/ , (2.45) 

due to the fact that f bac = —f cab in the basis (12. ip . In the equal coupling case K = /jwe can 
simplify the BPS equations (I2.43P and (I2.44p to the following combined equation 

- 3 ^{{HH<-,fl N ) a HH<) n . (2.46) 



'12 



Now we have seen the differences between the two different BPS systems in Sec. 12.11 and 
Sec. 12.21 The equation for the Higgs fields H (Eqs. (12.161) and (12.421) ). is common. The difference 
resides only in the flux equations (I2.19P and (I2.46p . which will eventually become the master 
equations after applying the moduli matrix method in Sec. HI We want to study the zero-modes 
of the BPS solutions to these equations in the subsequent sections. 
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2.2.1 The Abelian Chern-Simons solitons 



Let us recall the solitons in the minimal model by choosing N = 1 and Nf — 1, namely the 
Abelian CSH model [9]. The BPS equations are 



VH = 




(2.47) 



It is known that there are two kinds of solitons: i) the topological solitons and ii) the non- 
topological solitons. Let (H TP ,F^ P ) be a topological solution and (if NTP , F^ TP ) be a non- 
topological solution (several numerical solutions have been obtained in Ref. [9]). The topological 
solitons live in the Higgs vacuum while the non-topological solitons live in the unbroken (CS) 
vacuum. Therefore, the asymptotic behavior of the Higgs field is 

\H TP \-+v, |F NTP |-+0, as r^oo. (2.48) 

For topological reasons if TP must vanish at the center of the vortex but this is not the case for 
the non-topological soliton 

H TP , H NTP -»■ v' , as r -> , (2.49) 

where v' G (0, v) is a constant. The topological soliton has a quantized magnetic flux f !2.39|) 
whereas the flux of the non-topological one is a continuum 

$ TP = 2nk , $ NTP = 2n{k + a) , (2.50) 

with k G Z> and a > k + 2. For both topological and non-topological vortices, the integer k 
corresponds to the number of the zeros of the Higgs field H. 

It turns out that the dimension of the moduli space of the Abelian CS solitons is 

dim c A4 acs = k + a-1 . (2.51) 
Here a the integer-part of the real number a [9]. 



2.2.2 The U{2) Chern-Simons solitons 

To describe some of the characteristic features of the non- Abelian extension of the CS solitons, 
we will discuss the simple example of G' = SU(2) with Nf = 2 in this subsection. We also set 
k = fi so that the Abelian solutions automatically solve the non- Abelian equations 



VH = , 




(2.52) 
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Let us first figure out the vacuum structure of the model. There exist three different vacua: (2) 
the Higgs vacuum, (0) the unbroken (CS) vacuum and (1) intermediate vacua 



H w = r u , H w = r u , h® = r u . (2.53) 






The symmetry of the vacua is SU(2) c+ f in the Higgs vacuum, U(l) c x U(l) c+ { in the intermediate 
vacuum and finally U(2) c x SU(2){ in the unbroken vacuum. As already mentioned, the Higgs 
and unbroken vacua are unique but the intermediate vacuum possesses a continuum. By acting 
with the SU(2){ flavor symmetry, one can transform the above given configuration into a generic 
point 



H^=v 91 92 , 1^ + 1^ = 1. (2.54) 




Since the overall phase of <f>i, 02 is gauged, the vacuum manifold is SU(2) /U{1) ~ CP 1 . One can 
exchange the first and second rows in the matrix, which however is SU (2) c gauge equivalent to 
the given form. 

In the Higgs vacuum we can construct the topological vortex as follows 

H= ( HTF 0\ | diag(0,„) , as r^O, ^ 
y v J y diag (v, v ) , as r — > oo . 

Clearly, the color-flavor symmetry SU(2) c+ f respected by the vacuum is spontaneously broken 
to {7(l) c +f at the center of the vortex. This yields Nambu-Goldstone zero-modes which are the 
so-called orientational zero-modes in addition to the position (translational) zero-modes. Thus 
the minimal topological vortex carries the zero-modes [2T] 

Ml^ = C x CP 1 . (2.56) 

In the unbroken phase, we can construct the non-topological solitons as follows 

H = /P NTP 0W diagKO) , as r -> , 
y 0/ [ diag (0,0), as r -> oo . 

The symmetry of the vacuum U(2) c x SU(2){ is spontaneously broken to U(l) c x U(l) c+ { at 
the center of the soliton. The spontaneous breaking of the global symmetry leads to the 
Nambu-Goldstone modes SU(2)f/U(l) c+ f ~ CP 1 . For simplicity, let us choose the minimal non- 
topological soliton (corresponding to a — 2) which in the Abelian case has only translational 
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zero-modes. Clearly, with the above symmetry argument it indeed carries the same zero-modes 
as the topological soliton 

M^tp = CxCF i ) ^2.58) 

where the first C corresponds to the translation zero-mode. We will denote this the non-Abelian 
non-topological soliton. 

The intermediate vacua are interesting. We can put either topological or non-topological 
solitons. Let us first construct the topological soliton. A naive way is simply to embed the 
Abelian topological solution as 

/ H TP (A f diag (0,0), as r^O, 
H= -> { 6V ' ; ' (2.59) 

\ 0j I diag (v, 0) , as r — > oo . 

Since the symmetry at the vortex center is larger than that of the vacuum, there are no orien- 
tational zero-modes. However, this is not the most general solution. The generic solutions can 
be obtained by embedding not the Abelian solution with iVf = 1 but the Abelian solutions with 
Nf = 2, hence the semi-local Abelian vortex solution which have been studied in Ref. [13]. Let 
(Hf TP , i?f TP ) be the fields of the semi-local vortex solution in the Abelian theory. Then the 
generic solutions of the U(2) model can be written as follows 



rSTP ttSTP 

H=\ 1 2 \^{ 




as r — > 

(2.60) 

as r — > oo 



with v' being a constant which should be determined by the equations of motion. It is known that 
the semi- local vortex (Hf TP , i?| TP ) has a complex free parameter, say a G C. We can identify 
\a\ with the transverse size of the vortex while the phase is the relative U(l) phase between the 
first and second Higgs elements. When we send a to zero, the generic solution becomes the naive 
embedding solution (I2.59p . Thus we conclude that the moduli space of the topological solitons 
in the intermediate vacua is 

M s k l p = C x C . (2.61) 

The first C is the position and the second corresponds to the semi-local size modulus a. This 
embedding solution is easily extended to the U(N) case. As explained before, there are N — 1 
different kinds of intermediate vacua, see Eq. (12.351) . Obviously, in the vacua with m = 1 and 
N — m = 1 we can embed the Abelian semi-local vortex of the topological and non-topological 
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kind, respectively. In the rest vacua with m > 1 and N — m > 1 we can have non-Abelian 
semi-local vortices. So far, the semi-local vortices have been observed only in theories in which 
the number of flavors is greater than that of colors (iVf > N), in the case of G' = SU(N), see 
Refs. j22j|23]. Here we first show that the semi-local vortices also exist in the model with N{ = N. 
Furtheremore, if one likes, one can also consider systems with Nf < N, for instance Nt_= 2 and 
N = 3, and construct the semi- local vortex there by choosing an appropriate vacuum. |3 



There exists also a semi-local non-topological soliton in the Abelian theory with JVf = 2, see 

be a minimal solution (which corresponds to a = 2). Then we 



rSNTP^ 



Ref. [EJ]. Let (Hf , H\ 
can embed this into the non-Abelian model as follows 



H 



v 

SNTP 



#SNTP 




as r — > 



as r — > oo 



(2.62) 



The solution has also one complex free parameter so the moduli space is again 



mZ7 



C x C . 



(2.63) 



In Sec. HI we will study generic solutions and explicitly see the structure of moduli space by 
using the moduli matrix method. 

3 The index theorem 
3.1 Calculation of the index 

Let us first make a general calculation of the index. Later we will apply the results to the YMH 
and the CSH models considered above. 

Suppose the zero-modes 77 of a set of BPS equations are described by a Dirac-type equation 



Br] = , 




5 In G = SO, USp theories it is however different due to the aiready mentioned complicated vacuum structure, 
and the existence of the flat directions allows for semi-local vortices even for the minimal number of flavors in 
order to put the theory on the Higgs branch 2H1 HZ] • 
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where 77 takes a value in a vector-space generally described by matrices with rjf being an Nf x N 
matrix while ?] a is an iV x N matrix. The inner product on this vector space is defined as 

(l jV )= [ Tr\4 m + 477J , (3.2) 



while the adjoint operator of D is denoted by O^ and is defined as usual (4,077) = (DU, 77). 

In this section, we will calculate the index of the operator D of the Callias type [39J (i.e. the 
generalization of the Atiyah-Singer theorem to non-compact manifolds), which is formally 

X = dim (ker O) - dim (ker D 1 ") = dim(ker D f D) - dim(ker DD f ) , (3.3) 



where we have used that the kernel o: 
a slightly different index defined bjc 



DM is the same as D f D (DD"M. We will however calculate 



t { W) = tr ( , ^ ,) - tr ( f ,) , (3.4) 

which clearly gives back the index X in the limit M 2 — >- if the system possesses a mass gap, but 
in most cases it will turn out to be independent^ of M 2 , hence it will prove useful to evaluate the 
index in the limit M 2 — > 00. Subtleties arise in the case where the continuum part of the spectrum 
extends down to zero and one has to be more careful [5U] . This is because an unbounded resonance 
gives a certain contribution to the index. In such cases, we have to subtract the continuum part 
from X [9] . 

Let us now make some general statements for the form of the operators O, First of all, 
we decompose the Dirac operators as 

O = iD + K , D f = iD ] + K ] , (3.5) 

where D includes only the differential operators d, B and K, contain the remaining parts. 
This yields 

Dtp = A - Li, DD f = A - L 2 , (3.6) 



where we have defined 



A = -D ] D = -DD ] , (3.7) 
L x = -iD ] K - iK^D - K ] K , (3.8) 
L 2 = -iDK ] - iKD ] - KK ] . (3.9) 



6 We will adopt the notation where tr denotes trace over states as well as over the matrices while Tr denotes 
only a matrix trace. 

7 This can be proven for a compact manifold and the result proves to be the same if the fields are sufficiently 
well-behaved at spatial infinity. 
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We assume that A is diagonal and composed by Laplacians, while L± t 2 contain only linear differ 
ential operators. Expanding in terms of the large mass-squared parameter, we obtain 

M 2 



DtD + M 2 
M 2 



M 2 [P- 1 + P^LxP- 1 + P- X L X P- X L X P- X + ...] , 

M 2 \P- 1 + P- 1 L 2 P- 1 + P- 1 L 2 P- 1 L 2 P~ 1 + ..] , (3.10) 

where we have defined P = A + M 2 . All terms beyond the second order in Li will vanish in the 
large mass limit M 2 — > oo. 

Now we will pursue a technique used in Abelian systems [H] to simplify the calculation, which 
will demonstrate that the terms K^K and KK* do not contribute to the index. We can express 
the index in the following four terms as 



x = x 1 +x 2 +x 3 +x 4 , (3.11; 



with 



Xi = lim iM 2 ti [P" 1 {DK ] - K ] D - D ] K + KD ] ) P- 1 ] , (3.12) 

M 2 — >oo 

X 2 = lim Mhilp- 1 (KK^ - K^K) P" 1 ] , (3.13) 

X 3 = - lim M 2 ti[p-\D^K + K^D)P- 1 {D^K + K ] D)P- 1 ] , (3.14) 

X 4 = lim M 2 tr \p- x (DK^ + KD^)p- l {DK^ + X^p- 1 ] . (3.15) 

M 2 ->oo 

The trick now is that X 2 + X 3 + X 4 = which can be easily proven by making use of 

[p-\ D] = [P~\ D ] ] = 0, DP^D ] = D^P^D = -1 + M 2 p-\ (3.16) 
Hence, we can write the index simply as 

X = Xx. (3.17) 

This formula is considerably simpler than the starting point (13.111) . 

Let us calculate the index for a class of operators which includes our models 



?I3 Ld K i 

K2 iT>l Ah 



^Ld = B + i o A T , P a T dj = d - i [A T , o] . (3.18) 



where the gauge fields A, A of the gauge group G is a given background configuration. The Dirac 
operator acts on fields of the following form 



(3.19) 

5A T 
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where 5H is an N x Nf matrix-valued field in the fundamental representation of G whereas SA 
is an N x N matrix-valued field in the adjoint representation of G. K\ (i.e. an iVf x N matrix) 
and K2 (i.e. an N x TVf matrix) need not be fixed yet, so we leave them as variables until we will 
need the specific forms. It follows that with D = iD + K 

_ ( l Nt 8 \ _ /- o A T K x \ 

D ~ \ l N d) ' K ~ \ K 2 [A\o]) ' (3 ' 20) 

and its adjoint operator W = iD^ + K< is defined by 

, fl Nf d \ , f-oA T K\ \ , . 

D ] = f , K ] = \ 2 , 3.21 

V 1 N 5)' \ K\ [A\o]j 

where K\ 2 are the adjoint operators of K\p. The Laplacian operator reads 

A = -D^D = -DD^ = -3d . (3.22) 

Now by a simple calculation we can show that 

DK ] - K^D - D ] K + KD ] = v ; 



[dA T -dA T ,o] t 

independently of E 

We can now easily calculate the index by means of Eq. (I3.17P 

lim X=iM 2 I d 2 x tr (x IP" 1 (DK ] - K^D - D ] K + KD ] ) P^\x) 

M 2 —>-oo J 

j 1 Tr (OA - OA) j 



d 2 k M 2 



(3.23) 



(27r)^Ip + M 2 ) z 

= NfNu , (3.24) 

which turns out to be independent of M 2 and the Abelian winding number v is defined in 
Eqs. {235) and ( 12381) . for YMH and CSH theories, respectively. 

The result ( 13.24}) is quite general and can be applied to all the solitons under consideration 
in this paper. Namely, the index ( 13.241) is applicable to all the models considered in this paper: 
YMH and CSH theories with a general gauge group 

G = [17(1) x G']/Z no . 



8 Even though the index theorem calculation does not depend on the constants i^i,2, they are needed to derive 
the vanishing theorem, see below. 
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Next we want to show that the kernel of D' is indeed vanishing. This statement is usually 
denoted the vanishing theorem, and we will need to compute those case by case. The applicability 
of the vanishing theorem makes the index equal to the number of zero-modes of the system 
described by the operator D. 



3.2 Yang-Mills-Higgs theory 
3.2.1 The index for Yang-Mills vortices 

The index theorem in the case of G' = SU(N) YM theory was already made in Ref. [17] while 
it was made with a general group in Ref. [UJ, hence we will just briefly review the calculation 
here as we will need it later. 

In what follows, we set e = g, in order to make the expressions more compact. Our starting 
point is the BPS equations H 2 . 1 6 1) and (12. 19ft . Let {H, W,A,A} be a BPS solution. Then we 
consider small fluctuations around it as H — > H + SH, H* — > W + SH\ A — > A + SA and 
A — > A + SA. The fluctuations obey the following linearized BPS equations 

WfSH - SAH = , (3.25) 
iV.JA - WJA = e - {{SHH^) G + (H5H*) G ) , (3.26) 

where we have used that <SF 12 = 2i (P a SA — T> & SA) . We need to introduce a gauge fixing in order 
not to count non-physical degrees of freedom and it will prove convenient to choose [T7] 

2 



iV a SA + iV a SA= 6 -{- {5Hrf) G + (HSH^) G ) . (3.27) 



Combining Eqs. ( I3.26|) and (I3.27p . we can now write these equations in the form of Eq. (13. ip with 
the following Dirac operator 




r? = 1 , (3.28) 



^ - HT ) Dt= f ^ ^ 
:4(Ho) G iVl) ' \-2(Ho) G iT>l 

where H is just the complex conjugate of H while W is the Hermitian conjugate as usual. Since 
this operator is in the class of operators (I3.18p . we can immediately apply the result of the index 
fl3~24p . From Eqs. d2~T5|) and (l3T24|) . we find the index 

NtNk 

X — lim X = . (3.29) 

m 2 ->oo n 
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3.2.2 Vanishing theorem for Yang-Mills vortices 



In order to establish that the index does in fact correspond to the number of bosonic zero-modes 
of the BPS vortex solutions, we need to demonstrate that the dimension of the kernel of the 
adjoint operator is indeed zero 



dim (ker D f ) = . 
The question is if there exist normalizable zero-modes \ for W 

D f x = , X = 



(3.30) 




(3.31) 



where xj is an Nf-hy-N matrix and xj i s an N-by-N matrix and they obey the following 
equations 

,2 



AXf = ~jX,H , V aX , = -i2 (Xf#% 



(3.32) 



which can be combined into the following equation 



G • 



(3.33) 



In the Abelian case this is simply the Schrodinger-type equation 

^2 



-dd + —\H\ 



. 



(3.34) 



Since the Schrodinger potential is positive semi-definite, we can immediately conclude that x& — 
and this in turn tells us from the second equation of f 1 3 . 3 2 j) that also Xf needs to vanish on the 
entire plane. It follows that there are no normalizable massless modes in the Abelian system. 

For the non- Abelian case, the Schrodinger equation seems to be somewhat more complicated. 
There is however a nice trick to solve this problem, which is to write the square of the adjoint 
operator on a given state as a sum of squared terms. It was first introduced in Ref. p2] in the 
case of G' = SU(N), and was also applied to a general group in Ref. [U]. The starting point is 
to take the complex norm of IB^x = as 

= (D+x.^x) 
Tr 



c 



Tr 



l^fXf + X,H\ 2 + \iV aXa - 2 (Xf#%| 

PfXt I 2 + |^aXa| 2 + |Xa# I' + 4 | {XfH%\* + C 



(3.35) 
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where we have set e 2 = 4 and the cross term c has been rewritten into two total derivative 



terms by using the BPS equation T>H = 0: c = Tr id (xfH XL) ~ yx&Hx] j ■ Since we are 
interested in normalizable physical zero-modes, we require the fluctuations to vanish at spatial 
infinity (\z\ — > oo). Thus the total derivative terms vanish. We now want to show Xa = and 
Xf = from the conditions 

Z>fXf = 0, ^aXa = 0, Xa H = 0, ( Xi H^) G = 0. (3.36) 

The N-by-N^ matrix field H (background field) has the maximum rank TV almost everywhere 
except at the vortex centers. Thus we conclude Xa = from x&H = 0. For Xf we need a more 
elaborated argument except for the case of G = U (N) and Nf = N [17] . Let us explain Xf = 
for the case of G' = SU(N) and N { > First, we solve T>x\ = by a combination of the 
background field S and a holomophic matrix function Hq{z) by x\ = Hq(z)S. Then consider a 
holomorphic G' = SU(N) invariant operator (i.e. an iVf-by-iVf matrix) M = x\H = H (z)H (z). 
Note that x\ can be reconstructed from M, because H has the maximal rank in the broken 
(Higgs) vacuum. However, M vanishes at infinity because Xf should be normalizable, hence we 
require H Q (z)H (z) — > 0. The holomophic matrix H (z) = is a unique solution. Thus we 
conclude Xf — 0- F° r the generic simple group G', similar arguments were done in Ref. |47j . 

Now we have confirmed that has no normalizable zero-modes. So the index given in 
Eq. ( I3.29P indeed counts the correct number of the physical zero-modes. 



3.3 Chern-Simons-Higgs theory 
3.3.1 The index for Chern-Simons vortices 

Next we study the index of the BPS system in the non-Abelian CS theory. For simplicity, we 
set k = /i, hence we will be using Eqs. (I2.42p and (I2.46p . To the best of our knowledge, the 
index theorem for the CS vortices has been made only in the Abelian case in Refs. [HI [H]. With 
the formula (13.231) . we can easily calculate the index also in the non-Abelian case with a generic 
gauge group G'. The fluctuations of the BPS equations read to linear order 

iV5H-5AH = 0, (3.37) 

167T 2 / V 2 \ 

W 5A -iVSA= ( {5HH*) HH ] + {HH ] ) G 5HH ] 5HH j ) 

k \ 2 / G 

+ ^ r ((H5H^) G HH^ + {HH^) G H5H^ -—H5H^\ . (3.38) 
k \ 2 / G 



3 We thank K. Ohashi for this point. 
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We choose the gauge fixing condition as 

167T 2 / v 2 \ 

W SA + W5A = — ( {5HH ] ) G HH ] + {HH ] ) G 5HH ] 5HH ] ) 

k \ 2 / G 

Combining Eqs. (13. 38ft and H3 .391) . we obtain 



iVSA 



IfiTT 2 / V 2 \ 

'SHH^) G HH ] + (M% 5HH* - J 6HH ^ • ( 3 - 40 ) 



k \ /G 2 / 



Taking the transpose of the Eqs. (I3.37P and (13.401) . we can now rewrite the above linear 
equations as the Dirac-type equation (13.11) with the operator D defined as 



{HH T (Ho) g + H o (HH T ) G - ^Ho 




(3.41) 



G 

Although this Dirac operator takes a somewhat complicated form, it is indeed in the class of 
operators (I3.18p . Thus we can immediately apply the result of Sec. 13.11 also to this case and 
hence the index is given by Eq. (I3.24p . 

Let us explain the result for each type of vacuum. First for the topological solitons the result 
is very sound and the index is exactly 

NtNk 

1 = 1 = N { Nu = . (3.42) 

n 

Note that the index is always an integer and this result depends on the gauge group G' only 
through tiq. We will see that this index indeed counts the number of zero- modes. 

In the case of the non-topological solitons in the unbroken vacuum, the magnetic flux is not 
quantized, so 1 is not an integer 

~ NfNW + a) 

1 = -^-^ , 3.43 

n 

with k' G Z> and a G M + . However, since the index 1 counts the number of the zero-modes, it 
should give an integer. In order to obtain the correct result, we have to subtract contributions 
from the continuum due to unbounded resonances [50]. Thus the correct index is given by 

T= W + *), (3 . 44) 
n 

where a is the largest integer less than a (a — a < 1). k' is the number of topological solitons 
living "inside" the non-topological solitons and the number is limited by the amount of non- 
topological winding a. Note that this index still includes some non-physical zero-modes which 
change the magnetic flux We will return this point later. 
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When the topological and non-topological solitons coexist, the index is mixed as 

where the first term is related to the topological solitons while the second term is related to the 
non-topological solitons. Here we intend that the k vortices are topological in the (partially) 
broken vacuum, while the k' vortices are topological vortices living "inside" the non-topological 
vortices in the (partially) unbroken vacuum. 



3.3.2 Vanishing theorem for Chern-Simons vortices 

Next, we would like to make a vanishing theorem (I3.30p for in the CS case. The adjoint 
operator reads 



-2(Ho) G if>l 
The question is if there exist normalizable zero-modes \ 



(3.46) 



= , X = MA , (3.47) 

where xj is an N^-by-N matrix and xl is an N-by-N matrix. Unfortunately, we cannot use the 
same trick as we used in Eq. (I3.35P since the cross terms associated with the completion of the 
square do not vanish even by the use of the BPS equations. Therefore, we have to solve the 
problem in a more straightforward way. Let us write down the Dirac equation as follows 

V aX . = ^2( X tHi) G , (3.49) 

Taking the adjoint derivative T> & of the second equation, and using the BPS equations, allows us 
to combine the two equations as 

V a V aX , = ^^(x,HH^ G + (HH^ G x,-YX^HH^ . (3.50) 

This is a non-Abelian Schrodinger equation for \& with the background configuration A{ and 
H. This tells us that D f has a non-zero kernel if and only if the Schrodinger equation has a 
normalizable zero-eigenstate. 
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Unfortunately, up to now, no proofs have been given even in the simplest Abelian case but 
Xa — has been expected in general not to have normalizable zero-modes [HI [11]. For the non- 
Abelian cases, we have nothing to add beyond the Abelian case. We would in general expect 
that Xa = and Xf — holds and that the index (13. 24ft counts the complex dimension of the 
kernel of D and hence the number of bosonic zero-modes of the BPS configuration. A partial 
evidence for us to expect this is that we will explicitly realize the same number of the physical 
zero- modes for the topological solitons as is given by the index theorem result ( 13.421) . 

4 The moduli matrix method 

In Sec. El we have demonstrated the existence of a certain number of zero-modes for the BPS 
solutions by using a formal mathematical computation of the index of a Dirac operator describing 
the BPS configurations under consideration. Here in this section, on the contrary, we will realize 
the zero-modes in a more explicit manner. The so-called moduli matrix method [281 Ell E01 [2Q] 
is suitable for that purpose which was first developed to describe the topological solitons in the 
Higgs phase of supersymmetric YM theories without CS terms. 

Let us apply the moduli matrix for the vortex system [301 EDI EI] in order to solve the first 
of the BPS equations (I2.16P and (I2.42p (which is the same in both the YM and CS theories): 

H = S^Hoiz) , A = -iS^dS , S e G c , (4.1) 

where H (z) is an N-by-Nf complex matrix which is holomorphic in z. H (z) is called the moduli 
matrix [281 I3T1 l30l [20] and the decomposition of the gauge field is called the Karabali-Nair form 
[52]. S takes a value in the complexified gauge group G c . For instance for G = U(N) the 
complexified gauge symmetry is GL(N,C). Note that the one-form A is not pure gauge. The 
gauge symmetry and the flavor symmetry act on the new matrices as follows 

(S~\ H ) {U c S-\ H U f ) . (4.2) 

Furthermore, there is an equivalence relation, which is denoted ^/-equivalence and it acts as 

(S, H ) ~ V(z) (S, Ho) , V(z) e G c , (4.3) 

with all the elements in V(z) being any holomorphic function with respect to z. In order not to 
change the winding number (energy) of the solution, we should impose that the determinant of 
V(z) is non-zero. 
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That is all for the first equation. The moduli matrix Hq(z) has all the information about the 
topological solitons in the YM theories (311 [30] . It is conjectured in Ref. [26] that the same holds 
for the topological solitons in the (non-Abelian) CS theories. 

Since the classification of the moduli space for the non-Abelian vortices has been carried out 
in the case of YM theories, we will not repeat it here. In the case of CS theories, the results 
of Ref. [26] claim that the moduli spaces of the YM vortices apply also to the Higgs phase of 
the CS theories. However, in the case of CS theory - as already mentioned in Sec. 12.21 - there 
exist not only the Higgs phase but also the CS phase (unbroken phase) and the partially broken 
phases. In these cases, we cannot use the arguments of Ref. [26] and furthermore, as we will see 
later, the moduli matrix is not sufficient to describe all of the moduli parameters possessed by 
the vortices. 

After introducing the S field in Eq. (14. ip . it is natural to introduce a local gauge invariant 
field EU GUI EUJ 



tt = SS j , tt = utt . (4.4) 

Note that the complex matrix field S can be decomposed as S = sS with the Abelian part 
s G U(l) c ~ C* and the non-Abelian part S G G' c . In the right equation above, we have 
defined u = \s\ 2 and Cl = SS'. We will see that this Q will play a central role in the following 
sections. Note that Q transforms as Q — > VQV^ with respect to the ^/-transformation. Let us 
also introduce a "current" defined as [52] 

J = dfWT 1 . (4.5) 

This J transforms as a holomorphic connection under the ^-transformation 

J VJV" 1 + dVV' 1 . (4.6) 

In terms of J ', the magnetic field can be rewritten as 

F l2 = -2S- x dJS , (4.7) 

(with F® 2 = —2y/2Ndd\oguj and F 12 = —2S~ 1 B(dClQ~ 1 )S). The magnetic flux is expressed as 

$ = - f Tr[F 12 ] = 2Tr f 3 J = 2N [ ddlogu; , (4.8) 
Jc Jc Jc 

while the U(l) winding number is given by 

v = - [ ddlogu . (4.9) 



7T 



27 



Let us define a covariant derivative with the connection %J. For example, for an adjoint field 0, 
transforming as <fi — > V(j)V~ l , then we have 

V v <l> = d<f>-[J,<f>\. (4.10) 

T> v (j) transforms homogeneously under the ^-transformation. 



4.1 Yang-Mills-Higgs theory 



For the readers who are not familiar with the moduli matrix, let us review several examples (the 
details of this topic can be found in Ref. |20j). As the first example, let us review the Abelian 
vortex. The moduli matrix for k vortices reads 

k 

H (z) = a + a x z + a 2 z 2 + ■■■ + a k _ x z k ' x + z k = JJ(z - m) , (4.11) 

%=\ 



where we have fixed the coefficient of z k to unity by using the ^-equivalence relation. The set of 
complex parameters {ao, a±, . . . , ak-i} are the free parameters (the position moduli) of the BPS 
solutions. Hence, the moduli space of the k Abelian vortices is 



M k u(X) 



(4.12) 



The complex dimension is k, which is consistent with the index theorem result ( 13. 24ft . 

The second example we will review here is the k — 1 vortex in the G' = SU(N) model with 
Nf = N flavors. After fixing the ^-equivalence, one can obtain the following matrix 



Hn(z) 



\ 



(4.13) 



1 6tv_i 

Z-ZxJ 

The set of complex parameters {bi, ■ ■ ■ ,b N ^} corresponds to the inhomogeneous coordinate of 
the orientational zero-modes CP^ -1 while Zq is the position (translational) zero-mode. Hence, 
we obtain the corresponding moduli space 



■Mf/(Jv) 



C x CP 



7V-1 



(4.14) 



The complex dimension is N which is again consistent with the index theorem (13. 24ft . In this 
way, the moduli matrix provides us a realization of the physical zero-modes. 
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We would like to emphasize that the use of the moduli matrix is quite easy. One just needs to 
choose a holomorphic matrix which satisfies the conditions posed on the invariants with respect 
to G and it should be consistent with the chosen boundary conditions. For instance, any matrix 
Hq(z) whose determinant is a polynomial in z of degree k, generates k non-Abelian vortices in 
the case of G = SU(N). When we choose the form of the moduli matrix H (z) we should fix 
the ^/-equivalence in order to eliminate unphysical complex parameters. Once we have done this, 
one can easily read off the physical moduli parameters for the solutions. For a generic gauge 
group G, we have to take further conditions into account. We will not discuss those details here 
as they can be found in Refs. [HJ H7] . 

In order to solve completely the BPS equations, we now have to solve the second equation 
of the BPS equations (I2.16p . Unfortunately, this is not an easy task since no analytic solutions 
have been obtained even in the case of the simplest Abelian models. Solving the second equation 
is very important in the following two senses: 

1) We have to confirm the existence of the BPS solution for each given moduli matrix. If the 
solution does not exist, it means that the complex parameters in the moduli matrix cannot 
be moduli parameters. 

2) We have to confirm uniqueness of the solutions of the second equation. If we are not able 
to show the uniqueness, it means that there might exist additional zero-modes which are 
not included in the moduli matrix description. 

These problems have been resolved only in the case of the k vortices in the Abelian-Higgs model 
by Taubes [53]. But there is no available proof for the non-Abelian case^. 

The strongest statement in favor of our belief of existence of the solutions in the non-Abelian 
models is that several numerical solutions have been found explicitly in the literature. Hence, 
we are quite confident that the existence holds in all our models, although we currently do not 
have a rigorous mathematical proof. 

On the uniqueness problem, we have observed that the number of the physical complex 
parameters in the moduli matrix coincides with that of the index theorem (I3.24p in all known 
models [171 HZ]- Although this coincidence is in favor of our intuition, we still miss a direct 
relation between the uniqueness problem and the index theorem. We will partially solve the 
uniqueness problem of the non-Abelian vortices in YM and CS theories in Sec. O 

10 This is true on the infinite C-plane. See however Ref. [S3] for non-Abelian U(N) vortices on compact Riemann 
surfaces. 
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Let us rewrite the second equation of (I2.16P in terms of Q and the moduli matrix Hq(z) 

3d\ogu = -^Tr[v 2 l N -tt n- 1 ] , (4.15) 

BidCih- 1 ) = {^~ 1 ) G , , (4.16) 

with VLq = HqHI and Q = uCl. We have used the following relation 

(X) G , = Tr [Xt a ] t a = Tr [Xt a ] i a , t a = S^S' 1 = S^S' 1 , (4.17) 

where t a is an automorphism of t a . These equations are called the master equations for the YM 
vortices. The equal coupling case, g — e can for a generic group neatly be written as 

e 2 

didtttt- 1 ) = — ( v 2 i - n n~ 1 ) G . (4.18) 

Before closing this section, let us see the asymptotic behavior of Q. It is again sufficient to 
consider the minimal winding solution in the Abelian model. Hence we take H (z) = z and then 
u approaches flo/v 2 = \z\ 2 /v 2 . Plugging this into the master equation, one easily obtains the 
asymptotic solution 

uj = ^[l + qK (m e \z\)] , (4.19) 

v 

where K$ is the modified Bessel function of the second kind and q is an unknown constant 
parameter which can be determined numerically. 



4.2 Chern-Simons-Higgs theory 

As mentioned above, the first BPS equation is common for the YM models and the CS models. 
Hence, the moduli matrix method which solves the first BPS equation explained in Sec. 14.11 can 
be applied to the CS models without any modifications [26]. This is one of the significant features 
of the moduli matrix method, namely its application range is indeed quite large. 

Since we have already solved the first equation, the remaining task is to solve the second 
equation of (12. 42ft . In Ref. [26], the BPS equations were given for a generic gauge group G, 
but the master equations were provided only for G' = SU(N), SO(N), USp(2M). Here we will 
provide the master equations completely independent of the choice of G' 

^logw = -^Tr [ W 2 1 W - flpO- 1 ] Tr [O O" ! ] - ^Tr [(fi ^%, OoO" 1 ] , (4.20) 

Bidnh- 1 ) = -^Tr [v 2 i N - fioir 1 ] (noQ-% - ^ ((fW 1 )^ fW 1 )^ • (4.21) 
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The equal coupling case, /x = k can neatly be written for a generic group as 

1fi7T 2 

Bidnn- 1 ) = — ((v 2 i N - n n^) G n Q n- 1 ) G . (4.22) 

These are called the master equations for the CS theory. 

Concerning the master equation, we have the same problems as in the YM theories. Namely, 
the existence and uniqueness of the solutions. To the former problem, we are in the same situation 
as in the case of YM theories, the best argument in favor of the existence we have currently, is 
the various numerical solutions to the non-Abelian CS BPS equations The existence of the 
solutions however is an important future problem but it is beyond the scope of this paper. In 
the Abelian case, the existence however has been proved in the topological case [55] and in the 
non-topological case (for radially symmetric solutions) [56] . 

As for the uniqueness problem, we will partially solve it in the Sec. [5j But there is a big 
difference between the YM theories and the CS theories. As we will see in Sec. the solutions 
to the master equation have their own moduli parameters, when we choose a vacuum different 
from the Higgs vacuum at the boundary. 

4.2.1 Abelian Chern-Simons solitons 

We will now briefly describe the Abelian solutions [9] in the moduli matrix formalism with N{ — 1 
flavor. Starting with the vacuum configurations. There are two vacua: i) the broken vacuum and 
ii) the unbroken vacuum. The Higgs field is H = v (up to U{1) gauge symmetry) in the former 
case, so it is described by 

broken vacuum : H — 1 , s _1 = v . (4.23) 

Note that the V^-equivalence has been fixed by the choice H — 1. In the unbroken vacuum, 
H = 0, we choose the moduli matrix as follows 

unbroken vacuum : H = , s" 1 = 1 , (4.24) 

Here we have set s = 1 by using the ^-equivalence relation. At first glance, the U(l) gauge 
symmetry seems to be broken since s _1 transforms as — > gs^ 1 with g G U(l). But this 
transformation can be absorbed by an according ^-transformation, such that no symmetries are 
broken. 

Let us next consider k topological vortices. It is generated by the moduli matrix Hq = 
[\ i=1 (z — Zi). The number of moduli parameters is k and is in accord with the index f)3.24p . The 
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master equation determines uj as 



2 k 



ddlogu = ^Ll[\z- Zl \ 2 u~ l [ v 2 - ] [ \z - ztfuj' 1 ) . (4.25) 
K i=i \ i=i J 

We impose the boundary condition for uj in such a way that the Higgs field approaches to the 

Higgs phase \H\ —> v. With respect to u, this is equal to imposing the boundary condition 

k 

u ^ v~ 2 Y\_\z ~ Zi\ 2 , as \z\ ->• oo. (4-26) 

i=i 

Note that this boundary condition is unique for obtaining a regular solution. We would like 
to stress that all the moduli parameters of the topological vortices are included in the moduli 
matrix. Unfortunately, no analytic solutions to this equation even for k = 1 are known. The 
asymptotic solution of k = 1 however is 

Wfc=i =v- 2 \z\ 2 [l + qK (m KOO \z\)} , (4.27) 

where q is an unknown parameter which can be determined numerically. There is no difference 
in the asymptotic form of the topological solitons of the Abelian-Higgs model and the CS model, 
see Eqs. (I4.19P and ( I4.27p . This is because they are topological solitons in the broken vacuum. 
A tiny difference lies in the masses of the Higgs fields m e versus m KOO and in the numerical 
coefficients q versus q. 

Let us next consider the non-topological vortex with k Higgs zeros and magnetic flux $ = 
— j c F u = 2ir(k + a). The total magnetic flux is not necessarily an integer. In fact, a can be an 
arbitrary real number. We again take the moduli matrix H = Yl^=i( z ~ z i)- Then the master 
equation is also the same as Eq. (I4.25p . We solve this with a different boundary condition for uj 
consistent with H — > 0. It is determined by the total amount of magnetic flux 

= -2dd log u -> -2dd log \z\ 2{k+a) , (4.28) 

for \z\ — t- oo. Thus the desired boundary condition reads 

k 

uj ^ C^lzl^JJlz - Zi\ 2 , (4.29) 
i=i 

where C is a unique numerical constant. H asymptotically approaches the unbroken vacuum as 
H — > C/\z\ a . In this way, we have found k zero- modes {z{\ in the moduli matrix as in the case 
of topological vortices. However, in Sec. I5.2[ we will see that Eq. (14.291) is not the most generic 
boundary condition on u for the soliton solutions with a fixed magnetic flux $ = 2Ti(k + a). We 
will also see that there exist additional zero-modes which are not accounted for in the moduli 
matrix but on the other hand reside in uj. 
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4.2.2 U{2) Chern-Simons solitons 



Let us next explain the solitons in the U{2) CS theory (« = fi) with N{ = 2 in terms of the 
moduli matrix. 

We start by describing the vacua as we did in the Abelian case above. As we have explained 
in Sec. 12.2.21 there are three vacua in the case at hand: (2) the fully broken vacuum, (1) the 
partially broken vacuum and (0) the unbroken vacuum. The Higgs fields corresponding to these 
vacua are given in Eq. (12. 53 p . while the corresponding moduli matrices are given by 

r(2) _ „. { 1 U \ rr(l) _ „, / P\ rr(0) 





'do' 

Here a, (3 G C parametrize the vacuum moduli space i.e. CP 1 . We fix these parameters as 
(a, 0) = (1,0) by using ^-equivalence and flavor symmetry. For all the vacuum states, S is 
chosen to be the unit matrix. 

Let us next consider the solitons. Since we are working with k = fi, most solutions in the 
U(2) case can be obtained by the embedding of some Abelian solutions. For instance, the k = 1 
topological vortex can be obtained by 

{*-*> (A / ^ +l S,r° |2 b*(z-z )\ , 

H = , n=[ 1+161 1 M , 4.31 

V b ij \b(z-z y i + H 2 J 

where u s is the k = 1 solution of the Abelian master equation (14.251) with the boundary condition 
u s — > v~ 2 \z\ 2 . The U(l) winding number is not integer but half-integer v = - J £ ddlogui = 1/2. 
We have realized two complex moduli parameters z$ and b in the moduli matrix. The former is 
the position modulus and the latter is the orientational modulus parameter of CP 1 , see Ref. [32] 
for details. Hence, we have two moduli (N^Nv = 2) which is in accord with the index result 

Next we study solitons in the partially broken vacuum. We are interested in the semi-local 
vortex in this vacuum. The moduli matrix and fl are respectively given by 

U = (- °) , (4.32) 

where Zq again denotes the position. The complex parameter a is a combination of the transverse 
size \a\ and a relative phase Arg a, while u s is the solution to the master equation 

4vr 2 




dd log u s 



v 2 



z \ 2 + \a\ 2 



z \ 2 + \a\ 2 



(4.33) 
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This equation should be solved with the boundary condition 

w s J ' ' ' • 4.34 

v z 

_ 1 

Thus, the C/ (1) winding number is again a half integer v = i J c <9<91ogcj s 2 =1/2. Hence, we find 
again two moduli parameters in accord with the index result ( 13.24}) . In the Abelian-Higgs model, 
the zero-mode a of the semi-local vortex is a non-normalizable zero-mode. We suspect that a in 
the semi-local CS soliton is also non-normalizable. 

One can also consider the semi-local non-topological solitons in the partially unbroken vac- 
uum. We do not work it out in this paper, although it might be interesting. The reader who is 
interested in this topic can easily extend our moduli matrix method to this case by comparing 
our method and the results in Ref. 1131. 



Finally, we study the non-topological solitons in the fully unbroken vacuum. To this end, we 
choose the following moduli matrix and Q 

*-(::)■ n -(T:)' 

where u s satisfies the master equation (14.251) with k = 0. Hence u s behaves asymptotically like 

_ i 

~ \z\ 2a . The U(l) winding number is v = ^ J c ddlogu£ = a/2, which is a half of that of 
the Abelian case. The complex parameter b is a realization of the internal orientational moduli 
CP 1 ~ SU(2)/U(1) which we have explained in Sec. 12.2.21 Suppose that the Abelian solution 
has a = 2 (minimal choice) and hence provides one modulus parameter, then we have found two 
moduli parameters in accord with the index result ( 13.24}) . see however Sec. |5]for more comments 
on this. 



4.2.3 Higher rank gauge group: U(4) 

Let us next highlight the vortices in the higher rank gauge group. For concreteness, we take 
G = U(4) in this section with iVf — 4 flavors. We will explain the vortices not in the full Higgs 
phase but in the intermediate and unbroken vacua since the latter ones have not been studied in 
literature. 

Let us first consider the following intermediate vacuum 

H = v diag(l, 1,0,0) , (4.36) 

where U(4) c x SU(4)f is broken to U(2) c x S[U(2) i xU(2) c+i ]. We can put both the topological and 
non-topological solitons in this vacuum. For instance, the minimal winding (k = 1) topological 
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vortex can be generated by the moduli matrix 



CI = v 2 diag(wi, ui2, W3, U4) , Hq — diag(z, 1, 0, 0) . (4.37) 



The corresponding master equations are of the form 

4vrV 



ddlogwi = — — (1 - \zru] r )\z\ uT » ( 4 - 38 ) 
4?rV 

aaiogo; 2 = — -(l-wj^wj 1 , (4.39) 



aaiogw 3 = , (4.40) 

<9<91ogw 4 = . (4.41) 

The first equation is exactly the same as the master equation in the Abelian case. We solve it 
with the boundary condition 

uji ->■ \z\ 2 , \z\ ->■ 00 . (4.42) 

The second one is solved by U2 = 1 and the last two can be solved as 0^4 = const. We can set 
the constants to one using a ^-transformation. Thus we find 

D, = v~ 2 diag(wi, 1, 1, 1) ->■ v~ 2 diag(|z| 2 , 1, 1, 1) . (4.43) 

Let us decompose ft into the Abelian and non-Abelian parts 

ft = v~ 2 uf diag (u?,u)P,Ui*, u^A . (4.44) 

From this, we can read off the U(l) winding number and the magnetic flux 

v = ^- J dd\ogu\ = ^ , (4.45) 

$ = 2-kv x 4 = 2tt . (4.46) 



The moduli matrix in Eq. (I4.37P and the corresponding master equations (I4.38p - fl4.4ip can 
generate other vortices in the other vacua. If we solve Eqs. f)4.38l) - fl4.41l) with the boundary 
condition for uj\ 

\z\ 2+2a 

0J1 — > — — — , as \z\ — >• 00 , (4.47) 

with 1^2,3,4 = 1, corresponding still to k = 1, it however leads to the non-topological soliton in 
the intermediate vacuum 

H = v diag(0, 1,0,0) . (4.48) 
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Its U(l) winding number and magnetic flux are 



If- i 1 + a 
v — — / dd\o^uj{ = , (4.49) 

$ = 2ixu x 4 = 2tt(1 + a) . (4.50) 



Furthermore, if we solve Eqs. fl4.38p -f l4.4ip with the boundary condition for 



2/3 



loi -> |z| , w 2 -> — — , \z\ — > oo , (4.51) 

we have the topological vortex with the minimal winding number (A; = 1) and the non-topological 
soliton in the different intermediate vacuum 

H = v diag(l, 0, 0, 0) . (4.52) 

In this case, Q is decomposed as 

II / - —I —I 2 —I —I —I — i\ 

f2 = v~ 2 diag I u*u 2 4 ' w i 4a; 2 4 5 w i 4uj 2 4 5 w i 4c<; 2 4 ) • (4.53) 

From this, we can read off the U(l) winding number and the magnetic flux 



1 1 



v = ^-Jddlog = \ + ^, (4-54) 

$ = 2ixu x 4 = 2tt + 2nf3 . (4.55) 

Note that, although this form is very similar to the previous example, the two configurations are 
completely different. In the previous example there is a non-topological soliton with a topological 
soliton "inside" and hence flux 27r(l + a) whereas in the case here, there are both the topological 
vortex with the quantized magnetic flux 2n and the non-topological vortex with the flux 2irf3. 

The last possibility of boundary conditions for Eqs. fl4.38p - fl4.4ip is 

\z\ 2+2a Ul 2/3 
wi -> — - — , u 2 -> —jj- , \z\ ->■ 00 . (4.56) 

This leads to the two non-topological solitons in the unbroken vacuum 

H = diag(0, 0, 0, 0) . (4.57) 

The solution has U(l) winding number and magnetic flux as follows 

v = \fdd\og = + I , (4.58) 

$ = 2ixv x 4 = 2tt(1 + a) + 2vr/3 . (4.59) 
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4.2.4 Higher rank gauge groups: U(l) x 5*0(4) and U(l) x USp(A) 



The solitons in the models with G' = SO(N) and G' = USp(N) can also easily be worked out 
using the moduli matrix. To be concrete, let us consider the examples SO (A) and USp(4) in 
this section. G' = S0(A),USp(4:) can be dealt with on the same footing in the moduli matrix 
formalism. A tiny difference is the sign e = ±1 of the invariant tensor J given in Eq. (12. 221) . In 
what follows, we do not need to distinguish between G' = SO (A) and G' = USp{4) since all the 
equations turn out to be the same in these two cases. 

Let us consider the moduli matrix 

Q = v ~ 2 diag(wi, wi, u 2 , uj 2 ) , H = diag(-2 fe , z k , 1,1). (4.60) 
With this at hand, the master equation (I4.22p is simplified as 

'd<91ogwil2 \ _ 4vrV f(l- \z\ 2k io^) \z\ 2k u^% 



L-2 



. (4.61) 



Note that oj\ and U2 are decoupled since we imposed the special coupling relation K — jj,. 

The boundary conditions of determine the type of solitons and the vacua. If we take the 
boundary condition 

uji ->> \z\ 2k , as \z\ -> 00 , (4.62) 
with U2 = 1, we get k topological (coaxial) vortices in the full Higgs phase 

H = v diag(l, 1, 1, 1) . (4.63) 

The solution is decomposed as 

O -2 1/2 j- / 1/2 1/2 -1/2 -l/2\ /, c ,\ 

i2 = v w 1 diaglcj]^ , ^ , u x jLO-l ' 1 . (4.64) 
Thus we can read off the U(l) winding number and magnetic flux 

v = - [ dd\ogul /2 = - , (4.65) 

71 J 2 

$ = 2nu x 4 = ink . (4.66) 
If we choose a different boundary condition 

| z |2(fc+a) 

lui —7- — — , as |^| — >• 00 , (4.67) 
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with 0J2 — 1, it leads to the non-topological soliton in the intermediate vacuum 

H = v diag(0, 0, 1, 1) . (4.68) 
The £7(1) winding number and the magnetic flux can be read off as above 

v = - [ddlogco 1 / 2 ^^ , (4.69) 

TT J 2 

$ = 2nis x 4 = 4tt(£; + a) . (4.70) 
If we take a different boundary condition 



wi ->■ |2| 2fc , w 2 -> , as |z| 00 , (4.71) 

we get k topological (coaxial) vortices as well as the non-topological vortices in the intermediate 
vacuum 

H = v diag(l, 1,0,0) . (4.72) 

The solution is decomposed as 

n _2 1/2 1/2 ( 1/2 -1/2 1/2 -1/2 -1/2 1/2 -1/2 lM /, „ q n 

Thus we can read off the U(l) winding number and magnetic flux 

1 f _ ( 1/2 \j2\ k i (3 



$ = 2nu x 4 = 4irk + 4vr/3 . (4.75) 



The last possible choice of boundary condition is 

| 2 |2(A:+q) | z |2/3 

uji ->■ — — , w 2 -> — — , as |z| -)■ oo . (4.76) 

This generates non-topological solitons in the unbroken vacuum 

H = diag(0, 0, 0, 0) . (4.77) 
The solution has U(l) winding number and magnetic flux as follows 

v = i/aaio g ( Wl "W 2 )=^ + f . («8) 

$ = Itxv x 4 = A7r(k + a) + Anf3 . (4.79) 
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As in the case of U(A) CS theory, the non-topological solitons in U(l) x 5*0(4) and U(l) x 
U Sp{4) carry non-Abelian zero-modes associated with the spontaneous symmetry breaking. To 
see this, let us consider a minimal example with the embedding solution 

/ ntp mtp x diag(t/, v' 0, 0) at the origin 

# = diag(# NTP ,# NTP ,0,0) -»■ < . (4.80) 

I diag(0, 0, 0, 0) at the infinity 

The global symmetry respected at spatial infinity is S77(4)f which however is spontaneously 
broken down to C/(2) c+ f x SU(2){ by the soliton (the same is true for USp(4)). Thus the non- 
Abelian moduli space is given by the Grassmannian 

M olicn ^ lon = m - sm * Gr 4 , 2 . (4.81) 

Note that this is different from the orientational moduli SO(4)/U(2) (USp(4)/U(2)) of the min- 
imal topological soliton in the SO (4) (USp(A)) model. We have observed that the topological 
and non-topological solitons accidentally have the same orientational moduli SU(N)/U(N — 1) ~ 
CP N ~ 1 in the U(N) CS model. This is however in general not the case. 



5 On the uniqueness of the master equations 

In this section, we will try to solve the uniqueness problem of the master equations ( 14. 1 5ft - 
(I4.16P or fl4.20p - fl4.21l) . which we pointed out in Sec. HJ For that purpose, we will consider small 
fluctuations around the gauge invariant field Q 

n^n s + sn, (5.1) 

with Q s = S s Sl being a true solution corresponding to a given configuration Qq- If the master 
equation has the uniqueness property!"! <5f2 must be zero in the whole C-plane. So our goal is to 
confirm that 5Q = 0. 

5.1 Yang-Mills-Higgs theory 

In order not to introduce unessential complication to the following argument, let us consider the 
equal gauge coupling case e = g. The master equation in focus is Eq. f!4.18p . The fluctuations in 

11 Since, at the boundary \z\ — > oo, the solution O s is known in terms of the lump solution for the semi-local 
topological solitons (which is the statement that fl s is in the vacuum manifold at \z\ — > oo) and furthermore that 
this solution is an algebraic solution and unique, we know that even if there exist another solution f2 it has to 
obey the same boundary condition at \z\ — > oo and hence it follows that Sfl — > for \z\ —> oo. 
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Eq. (15. ip obey to the linear order 

dd ( j>+ [n^an^dt] = ^n; 1 (a 0^7% a , (5.2) 

where we have defined the N-bj-N gauge invariant matrix field 

= n^sn . (5.3) 

Note that since $ and deft transform as holomorphic adjoint fields with respect to the V- 
transformation 

f VtfV- 1 , 8^ VdtfV- 1 . (5.4) 
then the above equation can be further rewritten as 

(V v d<f>y = ^ (n^Qo^a , (5.5) 
where we have used the holomorphic covariant derivative (I4.10p and 

q; 1 (xq~ 1 ) g a = Tr [xn-H a ] n;H a n s = Ti [n^xi ] i a = (n; 1 x) G , (5.6) 
with t a = n~H a n s . 

We compare Eq. (15.51) with Eq. (I3.33|) by changing the variable from the adjoint field Xa to a 
gauge invariant field given by 

^ = S^S; 1 . (5.7) 
With respect to ip, Eq. ( I3.33|) is expressed as 

2 2 

(v v d^y = [ir^o^l t a = ~2 ( n ~ ln ^) G > ( 5 - 8 ) 

with t a = (Sl)~ 1 t a Sl. Thus, the fluctuation of the master equation and the zero-mode ip of 
the adjoint operator obey exactly the same equations (15. 5p and (15. 8p . 

Since we have already proven that ip = (% a = 0) in Sec. I3.2.2[ we can immediately conclude 
that = 0. Thus we have proven the uniqueness of the master equation. Then we conclude that 
all the zero-modes reside in the moduli matrix H and Q has no moduli parameters. 

Note that this uniqueness is local uniqueness since <fi is only a small fluctuation around the 
true solutions. We still do not provide any proof for the global uniqueness. This is one of the 
important future problems. See however the discussion in Sec. |6j 
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5.2 Chern-Simons-Higgs theory 



As in the YM case, we take equal couplings k — \i and consider a small fluctuation around the 
moduli matrix field Q — > Q s + 5Q which by plugging into the master equation (I4.22p yields to 
linear order 

where we have used the identity (12.451) before inserting the fluctuations and we have again used 
the gauge invariant field <fi = Q~ 1 5Q. As in the YM case we can use (15.61) to rewrite the above 
equation as 

(vM) ] = ^ (<^ <^o0) G ) G + (n: 1 ^ (n; 1 ^ - v 2 i N ) G ) G ) ■ (5.9) 

Next, we would like to compare Eq. (15. 9 p with Eq. (I3.50P which the zero- modes Xa of in 
Eq. (I3.46P obey, as in the case of the YM system. Using again the change of variables (15 .7p . we 
can write Eq. (13.501) on the form 

(2W) f = ^ (<^o <!^U^) G ) G + (V; 1 ^ (^o - ^) G ) G ) , (5-10) 

where (X) G = Tr [Xt a ~\ t a with t a = (S^)~ 1 t a S^ . Hence, the fluctuation of the master equa- 
tion and the zero-mode ip of the adjoint operator D' obey exactly the same equations, namely 
Eqs. O and f l5TT0|) . 

Although ip and <fi obey exactly the same equation, they are different fields. While <fi is 
a small fluctuation around Q s = S s Sl which is related to the physical field configurations H 
and Ai through Eq. (14. ip . ip contains the normalizable zero-modes of and has nothing to 
do with H and A t . This difference must be carefully dealt with, because <fi is not necessarily a 
bounded solution especially in the case of the non-topological solitons. What we should require 
the normalizability of, is not the small fluctuations of Q but of the physical fields Ai and H, 
namely &4j and 6H must be normalizable. 

To be concrete, let us first try to find additional moduli parameters of non-topological solitons 
in the Abelian CS case [UJ by taking advantage of the moduli matrix method 12 !. Given a moduli 
matrix Ho(z) = 1, we consider a fluctuation <fi = u~ 1 5u which asymptotically satisfies 

ddcf> = (2a;- 1 - „ 2 ) .-V - -^SrJ , (5-11) 



12 Here we do not fix any gauge unlike Ref. [TT] where the authors worked in the Coulomb gauge. Instead, we 
deal with the gauge invariant quantity, namely cu and its fluctuation. 
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where we have used uj ~ \z\ 2a /C. Furthermore, the variation in the magnetic flux density is 
given by 

\z\ 2a 

Fn + SF 12 = -2dd\ogu - 2dd<p -2dd\og^- - 2<9<90 , for \z\ oo . (5.12) 

In order not to change the total energy of the soliton, the second term must vanish. Therefore, 
the fluctuation will asymptotically approach the real part of the holomorphic function as 

/ 

<P^F (z)+F*(z) , F {z) = Y J ^ 1 - (5-13) 

i=i 

From Eq. (15. lip . / G N>i should be / < 2a. Note that we have suppressed the constant term 
(i = 0). This is because the asymptotic behavior of u is fixed as u — > \z\ 2a /C. For instance, if 
we consider F = a (a G C), we have to' = tu + 5tu — > tu(l + 2Re[a ]) ~ i+2Rc[a ] | z |2a_ gj nce we 
have chosen, however, the asymptotic behavior uj — > \z\ 2a /C in such a way that the configuration 
becomes regular, we have to choose do = 0. Finally, we should impose the normalizability 
condition on 5H and 5A. Fluctuations of H and A are given by 

5H = -FH , SA = -idF , F = s^Ss . (5.14) 

<fi and F are related by <fi = F + F*. Thus we can asymptotically identify F — > Fq{z) and 
in turn 8A — > is automatically ensured. Moreover, normalizability of 5H which behaves as 
5H — > —F \z\~ a requires that the power of F should be / = a — 1. Hence, we obtain the 
boundary condition of (ft and uj 

(5.15) 
(5.16) 

This supplies a — 1 additional complex parameters {a^} (z = 1, 2, . . . , a — 1) in uj as varieties 
of the boundary condition for oo. Note that this result is not consistent with the index theorem 
result X = a (iVf = N = n = 1) given in Eq. (13.441) . This mismatch was first observed in 
Ref. [9]. In this way, the index sometimes over-counts the zero-modes. So we should be careful 
when we count the physical zero-modes using the index theorem, i.e. we should correctly subtract 
unphysical zero-modes. This problem seems however only to arise in the non-topological cases. 

Therefore, we conclude that u has its own zero-modes in addition to those in the moduli 
matrix for the non-topological solitons. 
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This is in sharp contrast to the case of the topological CS solitons in which there are no 
degrees of freedom in u once the moduli matrix is given. This can be seen as follows. The 
fluctuation asymptotically satisfies 

dd(j> ~ m 2 K00 4> , (5.17) 

where we have used the asymptotic form of u ~ \z\ 2k /v 2 . Clearly, the right-hand side does not 
allow to be a harmonic function and there is no normalizable solution compatible with the fact 
that needs to be the real part of a holomorphic polynomial and hence we can conclude that 
must vanish. 

Once we fix the boundary condition of u, namely we set — > 0, we can utilize the similarity 
between and ip which both satisfy the boundary condition 0, ip — > at infinity. Since we have 
already shown that plausibly -0 = (% a = 0) in Sec. I3.3.2[ we expect in general that = 0, 
so that the master equation is unique. Note that this uniqueness is local uniqueness since is 
only a small fluctuation around the true solutions. As in the case of YM solitons, the global 
uniqueness problem is an important future problem. 

Let us next consider the non-topological non-Abelian CS solitons. As in the Abelian case, we 
consider a variation in the magnetic flux density 

5Tt[F 12 ] = -25 (<9<9Trlogft) = -2<9<9Tr[0] , (5.18) 

where = fi" 1 ^ is an N-by-N matrix field. The contribution to the total magnetic flux from 
must be zero, so we require the following boundary condition for 0: 

0^^ 1 Fo(z)fi s + F o t (*) , (5.19) 

where F (z) is an arbitrary iV-by-iV holomorphic matrix. With this boundary condition, we 
have to solve Eq. (15.91) . This holomorphic matrix F Q (z) supplies additional zero-modes to the 
non-topological solitons. In order to find a condition on Fo(z), let us consider the fluctuation of 
A and H which can be written as 

SH = -S^FSH , 8A = -iS^dFS , F = 8SS~ l . (5.20) 

Note that F is related to by 

= VL~ l FVL s + F f . (5.21) 

Comparing this with Eq. (I5.19p . we can identify F Q (z) as the asymptotic function of F, namely 
F(z,z) — > Fq(z) as \z\ — > oo. The fluctuations of the physical fields H and A must asymp- 
totically go to zero. The holomorphy of F Q (z) is indeed needed for SA to vanish asymptotically. 
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Furthermore, powers of the holomorphic functions in Fq(z) are determined by the normalizability 
condition on SH. 

Let us give some examples in order to understand better the situation. For simplicity, let us 
take G = U(2) and consider the non-top ological soliton in the unbroken vacuum H = diag(0, 0). 
To create the vortex, we take the moduli matrix and make the diagonal Ansatz for fl s 

H = diag(l, 0) , a = diag(wi, 1) . (5.22) 

We choose the asymptotic behavior u± ~ \z\ 2a /C as \z\ —> oo. Note that this moduli matrix 
does not completely fix the ^/-equivalence. The residual infinitesimal ^-transformation takes the 
form 

V = y ' . (5.23) 



,0 l + 5 2 (z 



This transforms fh as 



* a ^° : JV- (524) 

Now we are ready to look at the fluctuations at infinity. To this end, we write F Q (z) as 

F =[ fu{z) fu{z) ) . (5.25) 



J2l(z) f 2 2(z), 



Then the fluctuations are expressed by 



5il = FQ s + n s F^H fn + f ^ u ^ + f A , (5.26) 

\ + f*2 /22 + J 22 ) 

5H = -S-'FHo -+ - K /2fl1 °) . (5.27) 

V /21 0/ 

We impose the square integrability condition on the fluctuations 

Tr [SHSHi] -» Tr [ft^FofiJitf] = ^l/iil 2 + |/ 2 i| 2 • (5-28) 

Hence, we find that the most generic form of the fluctuations is given as 

&-i 

/n = 5>^, / 21 = 0. (5.29) 

i=i 

On the other hand, there are no constraints for / 12 and f 2 2, it seems that /i 2 and / 22 have an 

infinite number of additional zero-modes. However, they are unphysical because they can be 



44 



eliminated by the residual ^-transformation given in Eq. (15.24p . Hence, we have found at — 1 
additional moduli parameters in fu(z). Furthermore, we have found the orientational moduli 
parameter b given in Eq. ( I4.35p . So in total we have found a — 1 + 1 = a complex parameters 
for the non-topological solitons. Again, the result is not consistent with the index theorem result 
X = 2a given in Eq. ( I3.44p . At this stage, we are not certain if the index over-counts the number 
of physical zero-modes a by including unphysical modes or we simply did not exhaust all the 
possible zero-modes in the above moduli matrix calculation. We will leave this problem for a 
future work. 

Let us consider another configuration by taking the following moduli matrix 

H = diag(l, 1) , fi s = diag(o;i, u 2 ) ■ (5.30) 

We fix the boundary conditions as oj\ — >■ \z\ 2ai /Ci and 002 — > \z\ 2oi2 /C2, while we define a = 
ai + &2- Note that there are no residual ^-equivalence relations in this case. The fluctuations 
are written as (where F (z) still is given by Eq. (I5.25P ) 

5Q = FQ S + Q S F^ + + UJ * f A , (5.31) 

\wi/ 21 + w 2 /i2 w 2 (/ 2 2 + fh) J 

(-1/2, -1/2, \ 

W 2 J 21 W 2 J22/ 

As before, we impose the square integrability condition 

Tr [ffl"ffl"t] ^ Tr [^Fo^] = ^ (|/ n | 2 + |/ 12 | 2 ) + ^ (|/ 21 | 2 + |/ 22 | 2 ) . (5.33) 
This gives the upper bounds for f TJ (I, J — 1,2) 

i=\ 
a j — I 

fu = J2a? J h\ (J^J). (5.35) 

Here we did not count i = for /// because the constant part is fixed by the boundary condition 
uj — > \z\ 2ai jCj. Thus we found {a.\ — l) + a 1 + («2 — 1) + «2 = 2(di + d 2 — 1) moduli parameters 
in {fu}. Again, there is a mismatch with the index theorem result (I3.44p . see the comments 
above. 
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6 Conclusion and discussion 



In this paper we have studied the non-Abelian solitons in the M = 2 supersymmetric gauge 
theories in three dimensions; viz. the YM and CS gauge theories. We have used a common 
framework for investigation of the non-Abelian solitons which nicely works both for the YM 
theories and CS theories. We have found several new solitons; the non-Abelian non-topological 
solitons and the non-Abelian semi- local (topological/non-topological) solitons in the unbroken 
and the partially broken vacua which exist in the CS theories. The non-Abelian non-topological 
solitons have internal orientational moduli analogous to the non-Abelian topological solitons, 
which however is in general different from that of the topological solitons. Furthermore, to the 
best of our knowledge, the present paper is the first work which has found the semi-local solitons 
in the models with the number of flavors less than or equal to the number of colors. 

In addition to the minimal solution, we have also paid great attention to the multiple solitons 
and their zero-modes. To this end, we have made use of two methods, counting the number of 
zero-modes by the index theorem and realizing the explicit zero-modes by the moduli matrix 
formalism. Those for the YM solitons given in Sec. 13.21 and 14.11 are all known and we have 
simply reviewed them. What we learned from these review parts of the YM vortices is that the 
moduli matrix formalism is a rigorous tool for investigating topological solitons Indeed, the 
index theorem counting [T7j and the number of the zero-modes realized in the moduli matrix 
formalism [30] do in fact coincide. In order to comfirm that the moduli matrix formalism indeed 
provides us with all possible moduli parameters, we have to solve the long-standing uniqueness 
and existence problems of the master equations (14.151) and (14.161) . As to the uniqueness problem, 
we have made a progress. We have proven the local uniqueness of the master equation by finding 
a one-to-one correspondence between the small fluctuations of the master equations and the 
vanishing theorem. 

We have also tried to count the zero-modes by the index theorem for the CS theories, although 
we could not prove the vanishing theorem. Furthermore, in order to go beyond the parameter 
counting, we have also applied the moduli matrix formalism which has already been used to 
investigate the topological solitons in Refs. [26j [2T] . For the topological solitons, we have reached 
at the same level as the YM case, namely the coincidence of the index theorem result and the 
number of the zero-modes realized in the moduli matrix. Actually, we have found no differences 
both in the index theorem calcuation and the moduli matrix formalism for the topological solitons 



13 Needless to say, the moduli matrix formalism has been already established. Many results obtained using the 
moduli matrix formalism are summarized in the review 1201. 
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Fig. 1: The fluctuations of the master equation outside the radius R have to be infinitesimal, since the 
fluctuations have to vanish identically at infinity. Inside R they could become finite. However, we have 
shown in Sec. [5] that the infinitesimal variation has to vanish. 

between the YM and CS theories, except for some technical details. This fact strongly suggests 
that the moduli matrix formalism is also a rigorous tool for investigating the CS solitons. If this 
is the case, the moduli space of the topological CS solitons is completely the same as that of the 
YM theory which has been studied intensively, which was conjectured in Even the metric 
of the moduli spaces are coincident in the two theories (the dynamics is however different due to 
the CS term [21]), so there is nothing new about the moduli space of the topologial solitons. 

However, this is not the end of story. Proper features of the CS theory appear when we study 
the non-topological and semi-local solitons in the unbroken and the partially broken vacua. We 
have also made use of the moduli matrix formalism and showed it works well too. A remarkable 
contrast to the YM theory is that the gauge invariant field Vt has its own moduli parameters 
in addition those provided by the moduli matrix, namely, there are zero-modes in the master 
equation. We have found them in the varieties of the harmonic boundary condition for the 
master equation. Once H and the boundary condition of Q are properly fixed, we would expect 
in general that the master equation in the CS theory is unique at least locally. We have given 
several examples, G = U(2),U(4),U(1) x 50(4) and U(l) x USp(4) where the Abelian solitons 
can be embedded solutions in the special case of equal CS couplings. The generalization to the 
higher rank gauge groups is straightforward. 

We have attacked the long-standing problem of the uniqueness properties of the master equa- 
tions and we have proved completely independent of the gauge group that a small fluctuation of 
the gauge invariant field Q around the topological solitons is in one-to-one correspondence with 
the vanishing theorem of the index theorem calculation. 

There is one subtlety, but as we shall argue, physically not too worrying fact about this unique- 
ness calculation. In this discussion here, we will consider only the question of local uniqueness 
proved for the YM case and expected for the topological CS case. The reason for which we 
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only claim that the uniqueness is local is that the fluctuations of the master equation considered 
are just infinitesimal fluctuations. Hence, one would immediately ask, what if the fluctuation is 
finite? Is there then a possibility for the existence of another solution and hence a parameter 
governing the different solutions? Let us first recall that in the vacuum, the solution to the 
master equation is completely unique in all cases, i.e. it is the well-known lump solution (even if 
it is technically singular in the local vortex case). Hence at an infinite radius, the fluctuations 
must vanish identically. Now let us assume that there exist two solutions, VLa,b- The finite 
difference is necessarily zero at infinity. Consider now a finite radius R (see Fig. [1]) much larger 
than the typical scale of the vortex system /x. Let us assume that the finite difference could be 
non-zero within the radius R. However, from that radius till infinity the difference is expected 
to be very small, say of the order of l/(fiR) <C 1. Thus we expect the difference to be exactly an 
infinitesimal deviation from the true solution. However, our uniqueness calculation shows that 
if the vanishing theorem is satisfied, then there cannot be a non-zero fluctuation deviating from 
a true solution. The somewhat physical argument tells us that in this case, we can expect the 
difference Qb — to vanish. A more rigid demonstration of this difficult problem is indeed 
welcome and left for future works. 

We will now discuss what we were not able to do in this work. For the YM theory, the global 
uniqueness and the existence problems of the master equation are unsolved^. For the CS theory, 
1) the proof of the vanishing theorem, 2) the uniqueness and exisitence problems, 3) matching 
between the parameter counting of the index theorem and the moduli matrix formalism for the 
non-top ological and semi-local solitons, 4) generic solutions of the non-topological and semi-local 
solitons, especially with the intermediate relative orientational moduli and 5) numerical solutions 
for the non-topological and semi-local vortices for generic coupling k ^ \i. We left these problems 
for future works. 

A comment in store is about the equal gauge coupling choice which we have made in this 
paper. It causes no problem for the number of moduli in a solution. Our derivations of all the 
generic formulae done for equal gauge couplings can trivially be extended to the generic gauge 
coupling case. We have made this choice for simplicity and aesthetic beauty. However, the 
solution to a given master equation (CS or YM case) is clearly different in the case of different 
gauge couplings. Hence, for the realization of explicit solutions, it is certainly important to 
consider the generic gauge coupling case and verify the existence of the solutions for all (finite) 
values of the couplings and values of the moduli parameters. 



14 Clearly, the existence of the embedded Abelian solution exists provided the gauge couplings are equal. What 
we mean here is that the solution exists in general for all values of the couplings and moduli parameters. 
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Another interesting thing might be the realization of the low-energy effective theory on the 
non-top ological CS solitons. For the topological CS solitons it has been studied in Refs. [211 
|24"| [25] . Dynamics and interactions of non-Abelian solitons, especially non-topological solitons, 
which may depend on the orientations are also interesting open problems. 

Non-Abelian solitons in the CS theory with the YM kinetic term may have similar properties 
compared with the solitons in the current work. The D-brane construction of the moduli space 
is another interesting problem. Furthermore, the non-relativistic limit may be also interesting, 
especially due to the integrability found in the Abelian non-relativistic systems [57] . 
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